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The thermodynamics of the quantum Heisenberg antiferromagnet on a square lattice is revisited 
through a linearized spin- wave theory which is well defined at any finite temperature. We re-examine 
in details the temperature dependence of the free energy, the internal energy, the entropy and 
the specific heat. Most conclusions of the thermodynamics in previous studies can be reproduced 
in our linearized spin-wave theory. Specially, our calculation at low temperature T < J agrees 
quantitatively with the numerical Quantum Monte Carlo simulation and high temperature series 
expansions. 
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The two-dimensional (2D) quantum magnets is one of the main subjects in the modern condensed matter physics, as 
they are closely related to the unconventional superconductivity in the cuprates[l[ and the novel quantum spin- liquid 
states in frustrated lattices. The previous studies are mainly concentrated in the low-energy physics. However, some 
i— 1 1 interesting physics, such as the crossover behaviorQ in 2D quantum Heisenberg antiferromagnet (QHAFM) and the 
' qj . pseudogap physics in the cuprates,Q occurs at high energy scales and are still in controversial. One reason for these 
controversies is the lack of a reliable theoretical formalism well-defined in full temperature region. 

Recently we proposed a perturbation spin- wave theory for the 2D QHAFM on a square lattice. 0] Our theory is 
based upon Takahashi's modified spin-wave theory, where in order to study the finite temperature disordered state, 
he introduced a constraint of the loss of the local moment. [6| It is found that the 2D quantum antiferromagnet 
is well described at low temperatures by this modified spin-wave theory, yielding the results consistent to the other 
theories and numerical simulations. (5, [ToL In the modified spin- wave theory, a mean- field ansatz is introduced 
to decouple the spin- wave interactions. Although it gives a good description of the low-energy physics, it results in 
the failure of the theory at high tem per atures T ~ J, where an artificial phase transition occurs which obviously 
O . violates the Mermin- Wagner theorem. 

i ^ i ' In our perturbation spin- wave theory, the loss of the local moment is also included to fulfill the Mermin- Wagner 
theorem. The corrections from the spin-wave interactions are studied via a many-body perturbation method, which 
t-H ■ removes the mean-field divergence in the modified spin-wave theory. The temperature dependence of the uniform 
static susceptibility from our perturbation spin-wave theory agrees well with the quantum Monte Carlo simulations 
and high temperature series expansions. || It shows that our theory is well-defined at any finite temperature and is a 
good formalism for the 2D QHAFM. Moreover it is found that even the linearized spin-wave theory with the spin- wave 
interactions neglected can reproduce the main features of the uniform static susceptibility. Thus it encourages us to 
use the linearized spin-wave theory, on its reliability and simplicity, to study the thermodynamical properties of 2D 
QHAFM. 

In this paper, we present our calculations of the thermodynamical variables of a QHAFM on a square lattice, includ- 
ing the free energy, the internal energy, the entropy and the specific heat. Comparing with the numerical simulations, 
• • ■ it shows that the linearized spin- wave theory can also reproduce nearly all of the features in thermodynamics. At 
low temperature T < J, our calculations of the internal energy and the specific heat agree quantitatively with the 
quantum Monte Carlo simulations and high temperature series expansions. 
Vh ' The Hamiltonian of a QHAFM on a square lattice is given by 



(ij) 

where J > and (ij) denotes nearest-neighbor sites. Since the ground state is a Neel ordered state, we separate the 
lattice (with 27V sites) into two sublattice, where sublattice A (B) is denned for the up (down) spins in the Neel order 
state. 
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FIG. 1: The temperature dependence of the Lagrange multiplier A in our linearized spin- wave theory of QHAFM on a square 
lattice with spins S = 1/2, 1 and 3/2. 



Following Takahashi, we introduce a Lagrange constraint into the Hamiltonian, 6] 

H X = - Ml^f . (2) 

l£A,B 

The lagrange multipliers are assumed as fii = fi if I G A and —fx if / 6 B, and \i is re-defined as /i = JzS(X — 1) 
for convenience in our following discussion. This Lagrange Hamiltonian fulfills the Mermin- Wagner theorem for the 
loss of the local magnetization. It leads to a finite gap in the spin-wave energy spectrum and thus the spin-wave 
excitations behave as bosons with finite mass at finite temperatures. 

In our spin-wave theory, the spin operators are expressed by bosonic operators based upon the Dyson-Maleev 
representation, 

= a\, = (25 — a\ai)a,i,S 7 i = S - a\ai, (3) 

for the spins in sublattice A and 

S 7 =b j ,S+ = b}(2S-b]b j ),S* = -S + b]b j , (4) 

in sublattice B. With this spin representation, the Hamiltonian H = H s + H\ can be expressed as summation of 
quadratic and quantic terms in these bosonic operators. In our following discussion, we will neglect the quantic 
terms which represent the spin-wave interactions and investigate the thermodynamics of QHAFM within a linearized 
spin-wave theory. The quadratic linearized spin-wave terms in the Hamiltonian can be easily diagonalized through 
the Fourier and Bogoliubov transformations. The free energy per site is readily obtained as[5( 

T ^ . /_ . . /e k \\ 1 



£ In (2 sinh -Jz\S(l + 2S) 



N ^ V \2T)) 2 

k 

+ \jzS 2 , (5) 



where the spin-wave spectrum = JzS ^/X 2 — 7^ and 7k = cos ^k x cos hk y . z = 4 is the coordinate number. 

The parameter A is unknown and determined by the free energy. Minimizing the free energy with respect to A leads 
to a self-consistent equation for it,[f| 



s -5 + ^E«*(^)^j- < 6 > 



The numerical calculation of the self-consistent equation Eq.([6]) has been present in details in our previous paper. [B} 
Here we only cite the main results. The temperature dependence of the Lagrange multiplier A is shown in Fig. 
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FIG. 2: Thermodynamics of the QHAFM on a square lattice with different spin magnitudes, S = 1/2, 1 and 3/2. (a) is the 
free energy per site, (b) the entropy, (c) the internal energy and (d) the specific heat. 



[TJ At zero temperature the spin waves condensate to the momentum (0, 0) which results in a long-range magnetic 
order. [(| [l3| In this case, A = 1 + 0(l/N) and the spin-wave excitation is gapless. At low but finite temperatures, 
the long-range order is destroyed by the quantum fluctuations. [l2j Then A has an exponential behavior and decays to 

its zero-temperature value 1 as A = 1 + | ( j|#) 2 exp (— ^^j?™ ) . Here m = S — 0.19660 is the spontaneous ordered 
moment at zero temperature. The spin waves thus behaves as bosons with a finite gap A = \/A 2 — 1. It is this finite 
gap that leads to the vanishing of the long-range magnetic order at any finite temperature. At high temperatures, 
the self-consistent equation gives a weak dispersion for the spin waves with A = -j^g In (l + -|) . 

After the temperature variation of A is obtained, we have calculated the uniform static magnetic susceptibility in our 
previous paper. [5| Now let us turn our attention into the thermodynamical properties. Fig. [2] shows the temperature 
dependence of the free energy, the entropy, the internal energy and the specific heat of the QHAFM on a square 
lattice. The free energy per site is obtained from Eq.((S]). In low temperature region, 

F (T) = —zJzS (S + 0.157948) + cT 3 . 

where c is a ^-dependence constant. At high temperatures, it decreases linearly with temperature, 

1 



F = -Tin 



(l + S^l + S- 1 ) 5 ] +^JzS 2 . 



This high-temperature behavior has been obtained by Takahashi,Q where the last constant term does not exist. 
Based on the good agreement of our internal energy and specific heat with the numerical simulations as discussed 
below, we argue that the finite constant term in our theory comes from the linearized approximation. 
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The entropy shown in Fig. [2] (b) is is calculated from 

Sn = Jj 51 t( nk + l)^(nk + 1) - n k \nn k ] , (7) 



N 

k 



where n k — e e k ^ T _ 1 . At zero temperature, the entropy is zero, which is consistent to the second law of the thermo- 
dynamics. It shows a square law ~ T 2 as temperature increases. When temperature is high enough, it approaches a 
constant value 



S n = In 



:i + 5)(i + 5- 1 ) 5 



This value is different from In (25 +1). It is also due to the neglect of the spin- wave interactions in our theory. 

The internal energy is calculated from the free energy and entropy as E = F + TS n . At zero temperature for 
spin-1/2, our calculation gives E = —0.657948. It is quite close to the value —0.669494 which has included the 



spin- wave interaction corrections. [14j and —0.669437 in quantum Monte Carlo simulation. [151] It shows that the spin- 
wave interactions have a very weak contribution to the internal energy at low temperatures. When temperature 
increases, the free energy increases with an approximate T 3 -law and finally approaches a constant value at high 
enough temperature. This asymptotic behavior has also be shown in the high temperature series expansions. (l6j A 
slight difference is that there is constant value E = ^JzS 2 in our theory, but it is zero in high temperature series 
method. This difference stems from our linearized approximation. A detailed comparison of the internal energy from 
our linearized spin- wave theory with the quantum Monte Carlo simulation and high-temperature series study is shown 
in Fig. [3] (a). A good agreement at temperature T < JS(S + 1) shows that the spin- wave interactions has ignorable 
contribution to the internal energy at low temperatures. 

The specific heat can be calculated from C v = T^f- which leads to the following expression, 

1 v 4-(JzS) 2 ATf 
""WT^ sinh 2 (e k /2T) ' () 

Here ^ can be obtained by the derivative of both sides of the self-consistent equation Eq.©. The temperature 
dependence of the specific heat is shown in Fig. [5] (d). It vanishes at zero temperature and shows an approximate 
T 2 -law at low temperatures. It exhibits a broad peak at intermediate temperature T ~ JS(S + 1) and then decreases 
as bT~ 2 at high temperature. A more detailed analyses of our numerical data shows that b = 0.9075 for S = |. Both 
the low temperature T 2 behavior and T~ 2 decay at high temperatures are consistent to the high temperature series 
expansions. [17| The comparison with the numerical simulations is given in Fig. [3] (b) . It shows that our linearized 
spin-wave theory has grasped main features in the specific heat: a low temperature T 2 behavior, a smooth crossover 
broad peak at intermediate temperature and a high temperature T~ 2 decrease. A detailed comparison shows that at 
low temperature T < JS(S + 1), our linearized spin-wave theory agree quantitatively with the numerical simulations, 
which indicates the validity of the spin-wave formalism in studying the thermodynamics of the QHAFM. 



I. SUMMARY 



In summary, we have calculated the thermodynamical quantities of the QHAFM model on a square lattice from our 
linearized spin-wave theory. Our calculations include the free energy, the entropy, the internal energy and the specific 
heat. Although the neglect of the spin- wave interactions in our theory leads to some difference in magnitude compared 
to the numerical simulations, it has captured main features of the thermodynamics. The quantitative agreement of 
the internal energy and specific heat with the numerical simulations at low temperature T < JS(S + 1) indicates that 
the spin-wave theory is an appropriate formalism to study the 2D QHAFM model. 
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